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The Continuum Shell Model is an old but recently revived method that traverses the boundary 
between nuclear many-body structure and nuclear reactions. The method is based on the non- 
Hermitian energy-dependent effective Hamiltonian. The formalism, interpretation of solutions and 
practical implementation of calculations are discussed in detail. The results of the traditional shell 
model are fully reproduced for bound states; resonance parameters and cross section calculations are 
presented for decaying states. Particular attention is given to one- and two-nucleon reaction channels 
including sequential and direct two-body decay modes. New calculations of reaction cross sections 
and comparisons with experimental data for helium and oxygen isotope chains are presented. 
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Q ■ I. INTRODUCTION 

! 

New horizons are in sight in the field of nuclear physics as we move away from the line of nuclear stability. In the 
£S) ■ everyday experience we observe only a tiny fraction of the nuclear world, while the recent advances in observational 
' techniques reveal a hidden realm of extraordinary nuclear complexes. The term "exotic" is commonly used to highlight 
the unusual nature of newly discovered nuclear systems where structure and stability are governed by intricate interplay 
of quantum many-body structure and dynamics of nuclear reactions. Weakly bound nuclei and unstable resonances 
appear as important links in the chain of nuclear evolution in cosmos, their structure and properties are central for 
energy generation in stars and production of elements in the universe. Furthermore, the quantum objects of mesoscopic 
[ nature are common to many fields of science including but not limited to atoms and molecules, nanoscale condensed 
matter systems, atomic clusters, atoms in traps, and prototypes of quantum computers p|. In many applications to 
""^ ' mesoscopic systems one needs to understand and utilize the features of marginal stability and strong coupling between 
. discrete structure and continuum. 
1 1 ' Mean field along with corresponding shell structure is a starting point in the theoretical analysis of a quantum 
many-body system. The shell structure itself becomes exotic on the borderline of the continuum M S S 13 ■ Single- 
particle, pairwise and cluster excursions into the continuum become essential forming halo states and resulting in 
complex mixing of internal many-body states and continuum configurations. The problems of continuum shell model 
(CSM), or combining the description of reactions with the structure calculations, have been discussed almost since the 
dawn of the shell model, and summarized in the classical text However, only recently, fuelled by the discoveries of 
exotic systems, acute need in theoretical understanding, and growing computational capabilities, this subject received 
due attention, so that significant advances have been made during last years. In this work we concentrate on the 
version of the CSM 0, |g that is going back to the Feshbach projection formalism 0, 0, ITjl IT3 and, even much 
earlier, to the approach by Weisskopf and Wigner [13 and works in atomic physics by Rice [Tj| and Fano [la. Ilq. 
Alternative formulations, such as 0,ElEl|2(j, see a review [2l|, were also successfully developed recently. 

The specific attractive features of the approach discussed below are the natural unification of structure and reactions, 
full agreement with the results of the traditional shell model (SM) in the discrete spectrum, correct energy behavior 
of resonance widths and reaction cross sections near thresholds, self-consistent consideration of isotope chains, and 
exact unitarity of the scattering matrix. At this point we use the standard residual interactions adjusted in numerous 
applications of the conventional SM, although the problem of better interactions in the continuum remains open (the 
first steps in this direction were made in Ref. |22|). 

We organize the discussion here starting with the formal description of the CSM approach in the following section; 
this will be the formalism that we imply under the term CSM throughout this work. In Sec. [H] we also discuss 
mathematical details of the formulation, relation to observables and the unitarity of the scattering matrix. Sec. 11111 
is devoted to the detailed consideration of different parts of the interaction; we consider the limit of the conventional 
shell model one-body decay channels, and sequential and direct two-body channels. The realistic applications are 
shown and compared with experiment in Sec. IIVI 
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II. FORMULATION OF THE CONTINUUM SHELL MODEL 
A. Effective Hamiltonian 

In what follows we assume that the many-body Hilbert space is spanned by Slater determinants constructed from 
the single-particle (s.p.) orbitals \ j) in the mean field. Using the notations of secondary quantization we denote the 
s.p. creation and annihilation operators for discrete orbitals as &t and bj labelled by a combined discrete label j. 

For the continuum states we use s.p. operators b^(e) and bj(e) which are labelled with the discrete index j and the 
continuous s.p. energy variable e. If properly constructed for a mean-field potential, see below the discussion related 
to one-body reactions, these states are automatically orthogonal and form a complete set. In this work, however, 
we require only separate orthogonality of the bound states, [6^,6^,]+ = Sjf, and of the continuum states normalized 
according to [bj(e),b^,(e')} + = S(e — e') Sjj>. 

The full many-body space can be separated into two parts. The set of TV-particle bound states |1; N) bj N |0) 

forms the "internal" space, V, where the index 1 = {ji, J2 ■ ■ • j/v} labels the Slater determinant, the m-scheme 
representation in the SM terminology. The remaining space, Q, is the "external" continuum, i.e. many-body states 
which contain one or more continuum s.p. orbitals. The external many-body states |c, E) are labelled by total 
continuum energy E and the set of asymptotic variables c that defines a reaction channel (it includes the characteristics 
of the residual nucleus). The channel variable c is discrete only in the case of the one-body decay where energy 
conservation fully determines energies of the two decay products. In general this variable is continuous, containing 
relative energies or momenta of decay products needed for full specification of the final state. In the description of 
the formalism below we use notations and S cc > which in the case of the continuous channel variables should be 
interpreted as J dc and 5(c—c'), respectively. In Sec. IIII Dl where channel c in the two-body decay implies the presence 
of the relative energy variable, the sum over channels is explicitly given in terms of an integral. By construction, the 
many-body states in V- and Q- spaces are mutually orthogonal and normalized as 

(1|2>=<5 12 , { C] E\c';E')=5 cc ,8{E-E'). (1) 

Within the total space V + Q, we solve the stationary Schrodinger equation 

H\a\E) = E\a;E), (2) 

where the full wave function \a; E) is in general a superposition of internal states |1) and external states |c; E') 1 

\a-E)=J2^i(E)\l)+J2 I dE'a c (E'-E)\c-E'). (3) 

1 c J 

The Hamiltonian H has parts acting within and across V and Q spaces, H = Hpp + HpQ + Hqq. The external states 
can be eliminated by introducing a propagator that acts exclusively within Q space, 

Gqq(E) = 1 (4) 

E — Hqq + iQ 

where the infinitesimal imaginary displacement selects the appropriate boundary conditions for the scattering problem. 
Next we assume that the channel labels c correspond to the eigenchannels (23] in the Q space, Hqq\c; E) — E\c; E). 
Then the Schrodinger equation (J2J "projected" into the subspace P becomes 
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a 2 = 0. (5) 



The amplitude of the continuum admixture in the full wave function J2Jl is 



ac(Sli?) - E-E> + l ' (6) 

where we introduced notations for the H-pQ coupling amplitude 

Al(E',E) = (l\H-E\c;E'). (7) 

This amplitude depends on the continuum variable of energy E' and running energy E. As follows from the definition 
Q, there is no explicit independence when internal and external spaces are orthogonal. However, the important 
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independence remains; the kinematic factors included in the definition of the channel states \c;E') ensure that the 
phase space shrinks to zero and the channel c becomes closed below threshold energy E c characteristic for a given 
channel. Whence, the amplitudes A\ vanish at E' < E c . 

The set of equations (jSJ for coefficients ot\ looks as an eigenvalue problem with the effective Hamiltonian matrix 7i 
in the intrinsic space defined as 

<ii«(25)i2> = am +E / dE ' ACl{E ' E E) f+^ E) - (8) 

The integral in Eq. JSJ) can be further decomposed into its Hermitian part (principal value), A(E), and the anti- 
Hermitian part, -\i/2)W{E), 

f , AUE',E)A C *(E'E) ^ f . A?(E' E)Af(E', E) ^ s Ank , s 

T.J dE E-E' + iO =J2 P - v -J dE E _ E ) ~™ E A\{E)Af{E), (9) 

c c " c(open) 

where A\{E) = Af (E, E). Thus, the effective Hamiltonian for the V space (|10() takes form 

H(E) = H vv + A(E) - % - W(E) , (10) 

The resulting dynamics generated by the effective Hamiltonian (|10|) contains a usual P-space contribution that 
we identify here with the traditional SM corrected by the virtual "off-shell" excitations into the continuum via the 
self-energy term, 

(1\A(E)\2) = £ P.v. / dE> ^(E',E)AnE^E) ^ (n) 

C 

and supplemented by the anti-Hermitian term, 

{l\W{E)\2) = 2tt At(E)A-*(E), (12) 

c (open) 

that represents the irreversible departure into Q-space, i.e. decays; see Fig.^ The term W that comes from the poles 
in integration JSJ), has a factorized form, which is shown below to relate to conservation of probability and unitarity 
of the scattering matrix. The amplitudes A C (E) here represent on-shell processes that depend only on one energy 
parameter and correspond to real decays with energy conservation, E' = E. These poles in integration appear only 
when running energy E is above the decay threshold E c for a channel c. The channels where decays are allowed are 
referred to as open. At E — > E c + the amplitudes A\ vanish due to the kinematic factors implicitly included into 
their definition. 

c ,< c 
aV"-.a c 2 aV 
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FIG. 1: (Color online) Diagrammatic equation for the full Hamiltonian corresponding to the dynamics in the V space, Eq. 
JTUJ. Parts H° and V on the fi gure indicate the one- (s.p. energies) and two-body parts of the internal Hamiltonian H-pv- 



B. Scattering matrix and unitarity 

The "outside" view from the reaction side of the problem is equally important. In accordance with general scattering 
theory, the transition matrix, 

T^m^Tm^E^m)^^ (13) 



describes the process that starts in the entrance channel c' with amplitude A\ originated from the interaction H-pQ, 
continues through internal propagation within the P-space driven by the non-Hermitian energy-dependent effective 



4 



Hamiltonian (|10|) (with all excursions into Q space included), and ends by exit to the channel c described by the 
amplitude A1*(E). The scattering matrix can be written as 

S CC \E) = exp(i£ c ) {S cc ' - 2iriT cc '(E)} exp(i^). (14) 

The additional phase shifts £, C (E) describe the potential scattering or a contribution of remote resonances outside of 
the valence space of the model. 

The factorized nature of the non-Hermitian contribution to the effective Hamiltonian is the key for conserving the 
unitarity of the S matrix [24[ • This can be demonstrated by considering the propagator for the effective Hamiltonian 

0(E) = ^ (15) 
generated from the unperturbed propagator for the full Hermitian part, 

G ^ = E-H V 1-A(EY ™ 

With W = 2?rAAt, where A represents a channel matrix (a set of columns of vectors A \ for each channel c), we 
iterate the Dyson equation 

g = G-(i/2)Gwg (17) 

and, due to the factorized form of W, come to 

Q = G — inGA . 1 Atg, (18) 

1 + ittAIGA 

that is called the Woodbury equation in mathematical literature. 

The transition matrix, T — A^QA, can then be written with the aid of the matrix R = A^GA that is analogous 
to the i?-matrix of standard reaction theory. The unitarity of S'-matrix follows directly from these equations, see also 

m, 

S^ 1 —^-. (19) 



1 + inR 1 + inR 



C. Energy dependence and resonances 

The effective Hamiltonian of Eq. i|l(J|) is energy-dependent: at each scattering energy E its running eigenvalues are 
complex numbers £ a {E). This highlights the structure in Eq. J3J that the eigenstate is a superposition of internal 
states and asymptotic decay states that have right energy. The relatively small, and numerically tractable, dimension 
of basis states, in exchange for non-Hermicity and energy dependence, is a noteworthy advantage of this method as 
compared to direct discretization of continuum used in other approaches. 

The eigenvalue problem involving a complex matrix of a general form requires finding two sets of adjoint eigenvectors: 
left, \a), and right, \a). They satisfy 

H\a)=£ a \a) and (a\H = £„{a\. (20) 

The left and right eigenstates correspond to time reversed motions; they no longer have to coincide because the T- 
invariance in the internal space is broken by irreversible decays. The global symmetry with respect to the direction of 
time is however maintained by the full Hamiltonian that includes the products of reactions. As a result, the left and 
right eigenstates have the wave functions interrelated by complex conjugation which is the time inversion operation. 
The Hermitian conjugation of the Hamiltonian switches the roles of left and right; the same effect can be reproduced 
by selecting an advanced propagator boundary condition in I0J discussed earlier. The biorthogonality relation is given 
by (a |/3) = 8 a p; similarly, the expectation value of an operator X is (a|X|/?). These properties, the formalism of CSM 
and its interpretation become transparent in the one-body problem discussed below in Sec lIII^I for further notes on 
this topic we refer to [2|J. 

There are various interpretations of eigenvalues £ a (E) of Eq. © . Only for bound states below all decay thresholds 
a condition £ a {E) = E is satisfied since W = 0. With non-zero W, the eigenvalues of the effective Hamiltonian are in 
general complex, 



£ a (E) = E a (E) - l -T a {E), 



(21) 
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describing the quasistationary states. These resonances and their widths r Q > satisfy the Bell-Steinberger relation 

(a\W\a) = T a , (22) 

where the left hand side can be expressed through the amplitudes A\ transformed to the biorthogonal basis of 
quasistationary states |a). 

The continuation of the original problem to the lower part of the complex energy plane, E — > £ = E — (i/2)T allows 
the condition 

£ a {£) = £. (23) 

The complex energy roots here can be identified with the poles in the scattering matrix 11411 , and the states correspond 
to many-body resonant Siegert states [53 since by construction the eigenstate Q is a regular function with outgoing 
asymptotics. One can also use the Breit-Wigner approach psj and identify resonances differently, with a condition 

Re[£ a (E)} =E, T a = -21m[£(E)\. (24) 

In the limit of small imaginary part (narrow resonances), both definitions are equivalent. However in general the 
difficulty in parameterizing the resonance width and centroid energy is related to the non-exponential character of 
decay caused by the energy dependence of the Hamiltonian parameters. Wide resonances cover broad regions of energy 
and therefore are particularly affected by this dependence. This leads to the non-generic and asymmetric shape of 
the resonance cross section that makes standard Breit-Wigner or Gaussian parameterizations inappropriate. 

With either definition, Eqs. I|23|) or l|24l) for resonant states are complicated sets of nonlinear equations. In some 
cases @ the condition (|23J) may lead to unphysical solutions. For the realistic calculations shown in Sec. II VI we select 
the Breigt-Wigner definition (|24(l and implement an iterative approach starting from the energy determined by the 
conventional SM with W — 0. Although it is convenient to express the solutions in terms of resonant states, the 
parameters are definition-dependent and become misleading for broad states or in the case of overlapping resonances 
when interference is important. In these cases one should turn to the observable scattering cross section determined 
by the S'-matrix of Eq. 114J1 . The computation of the scattering cross section is a problem of matrix inversion 
which is linear in accordance with physical principles, but it has to be done at each energy which can make this task 
numerically unstable for narrow resonances. The cross sections can be calculated using the i?-matrix of Eq. 119|) and 
the Woodbury equation in which case complex arithmetics can be avoided. We will see complementary pictures that 
can be obtained by calculating the cross section and via resonance parameters coming from the diagonalization of the 
effective Hamiltonian in the example shown in Sec. II VI 



III. FROM HAMILTONIAN TO DYNAMICS 



The derivation above is based on the decomposition of the full Hamiltonian into H — H-pv + Hvq +Hqq. Another 
useful classification traditional to the SM is by the types of many-body processes it can generate. At this stage we 
restrict the interactions by one- and two-body and limit the space Q by the states with only one or two nucleons in 
continuum. The typical shell model limitations by few valence shells are imposed on the intrinsic space V '. In this 
framework we define the full Hamiltonian and discuss processes associated with each of the terms. 

Our discussion here, however, does not touch the lack of knowled ge o f the effective interaction. Although sophisti- 
cated methods of deriving the effective interactions were suggested [23, [2^ , the best results and the most predictive 
power in the conventional SM come from the phenomenological interactions, such as USD [30| , fitted to experimental 
data. The situation becomes increasingly more complicated when effective interactions involving continuum are to 
be used j2(J . In this work we identify the internal interaction H-p-p , together with the Hermitian self-energy term A 
included, with a shell model Hamiltonian of the form 

H vv + A = 'A 1 ', + \ E nhh)hH)b]b]bj*K- (25) 

3 

The parameters, s.p. energy levels and antisymmetrized two-body matrix elements, are known from fits to experimen- 
tal data, such as [3(j, see also review 0, and available from interaction libraries, such as j^. The energy dependence 
of these parameters that comes from A(E) is ignored here since it has not been included in the fitting process. As 
demonstrated below this dependence is weak and smooth function of energy. Defining the internal Hamiltonian in 
this way we guarantee that below thresholds the CSM provides the results identical to the well established SM with 
effective interactions. Above thresholds, the shell model interactions were fitted to experimental data using i?-matrix 
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analysis that identifies interactions (|25|l in the same way, Eq. I|19|) . The consistent readjustment of shell model in- 
teraction parameters taking into account energy dependence of Hermitian part is beyond the scope of this paper but 
remains a subject of future work. 

In this work we assume that the Hamiltonian describing the motion of nucleons in the continuum is purely single- 
particle 

3 

The asymptotic one- and two-particle continuum states considered in this work are antisymmetrized products of 
internal eigenstates a of the residual nucleus and the wave function(s) of particle(s) in the continuum. The states in 
one-body channels, 

\c,E) =b](e j )\a;N-l), E = E a + e j} (27) 

are labelled by energy E and the discrete channel index c that combines the s.p. quantum numbers j and characteristics 
a for the eigenstate of an (N — l)-particle daughter system. 

The assumption l|2tj|) allows one to express similarly the two-nucleon channel states, 

\c;E)=b\(e)bt,(e')\a;N-2), (28) 

characterized by the total energy E = E a + ej + combined of the daughter binding energy E a and energies of 
emitted nucleons in s.p. states j and j'. Here the channel c contains in addition information on continuous relative 
energy distribution between emitted particles. Within this work we do not consider cases when both particles in the 
continuum are charged thus the assumptions in Eqs. (|26[) and (|28|l are sufficient. The generalization for bound states 
of two particles in the continuum is relegated to the next stage. 



A. Single-particle decay and resonances 

Now we need to define the interaction Hvq responsible for coupling of intrinsic space with the continuum. We start 
with the one-body part and associate it with some potential V assuming this potential to be spherically symmetric. 

We first treat a pure s.p. problem of a particle moving in the mean field or independent particle shell model. 
The generalization covering all s.p. channels in many-body cases is straightforward and discussed in what follows. 
Of course, this subject is extensively covered by textbooks [32L l33l l34|. The purpose of the formulation presented 
here is to emphasize the conceptual identity between the full CSM and its trivialized version represented by a single 
particle in a potential. The notions and definitions of resonances, scattering matrix and its poles, time reversal and 
non-Hcrmiticity already appear in this simplest case. This section also highlights some technical details used later for 
the s.p. part in the full CSM including the generic threshold behavior of the decay amplitudes. 

In the coordinate representation the Schrodinger equation for the radial part of the s.p. wave function, 

(r\b]\0) = [Y lX ] j ^-, (29) 
where Yi and \ represent the angular and spin parts coupled to total angular momentum j, is 

} Uj (r) = k 2 Uj (r), (30) 

where k 2 = 2[ie, [i is the reduced mass, z and Z are charges of the particle and of the residual nucleus, respectively. 
The spin-orbit part can be included here assuming that the potential V{r) depends on / and the spin orientation 
which in our notations are hidden in the s.p. index j. 

For Eq. I|3(J|) with V(r) — 0, we have the regular, Fi(kr), Fi(0) — 0, and the irregular, Gi(kr), solutions as Coulomb 
wave functions with the charge parameter r\ = fiaZz/k. For a neutral particle, z = 0, the regular and irregular 
solutions can be expressed in terms of spherical Bessel and Neumann functions, 

Fi(kr) — krji(kr) , Gi(kr) — ~krni{kr). (31) 

The two independent solutions are related by the Wronskian, 



1(1+1) 



dr 2 



2[i 



V(r) + e 2 — 
r 



(32) 
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Thus, the Q space states are energy-normalized regular solutions, 



{r\j;e) = (r\b](e)\0) = [Ya] j M^. 



(33) 



Following the definition in Eq. (JJJ), the s.p. decay amplitude is 



°j( c . 



j 1 



(j\H-pQ ~ e|i; £j) 




drFi{kjr) [V(r) 



e] u,-(r). 



(34) 



A positive energy internal state decays with the width 7j = 27ra 2 determined by the above equation; under this 

choice of phases, the decay amplitudes are real. The s.p. amplitude has only one index j. 

This result for the decay width can be reproduced through the equivalent consideration of the on-shell scattering 
process. Let us introduce incoming and outgoing (Coulomb) waves Of{r) = Gi(r)±iFi(r). Consider a resonant state 
Uj] since through the rest of this subsection we concentrate on a state with a given set of s.p. quantum numbers j 
we will omit this subscript in notations, the orbital momentum subscript I which is a part of the combined index j is 
also omitted. The state u can be normalized as a discrete state when the decaying component that obeys the Siegert 
|27| outgoing wave boundary condition, 



lim u(r) =AfO + (kr), 



(35) 



is neglected. Then the outgoing flux normalized by velocity determines the decay width 

7 = 2™ 2 = -|AA| 2 . 

It follows from here that the asymptotics of the decaying states are given by the decay amplitude, 



(36) 



lim u(r) = —\j —j^- a(e)0 + {kr) 



(37) 



where we selected a phase to be consistent with the previous definition. Using the Wronskian relations, the outgoing 
part can be extracted from the wave function Uj leading to 



o(e) = - 



1 f dF 
u — — b — 

litiik \ dr dr 



(38) 



This equation is identical to Eq. I|34|l since the Schrodinger equation (|30|l that must be used to determine the outgoing 
component guarantees that 



d ( dF dw\ 

T (u^-F-)=-2»FV(rMr) 



(39) 



where F is any of the Coulomb wave functions. 

The eigenstate wave function in the asymptotics can be expressed via the s.p. scattering phase shift, Uj(e) 
cos(Sj)Fi + s'm(5j)Gi. The related ^-matrix then can be found as 



S = exp(2iS) = 



i4-o--o-±i 



u-fO+ ~0+-fu 

dr ar 



(40) 



that is consistent with definitions (|13fl and l|14|) . The poles of the scattering matrix correspond to the condition of 
the regular wave function with the outgoing wave in asymptotics. Just as in the general case in Sec. [H] this can 
not be satisfied at real energy, while if the problem is taken into a complex energy plane, k — > n = k — ik and 
e — * e = e — ij/2, the discrete set of solutions emerges. Thus, the resonance energy e and the width 7 can be defined 
as real and imaginary parts, respectively, of the complex energy which is the pole of the scattering matrix. Consistent 
with the general theory, the time reversed problem is physically equivalent; the boundary condition then is that the 
wave function be regular at the origin and represent an incoming wave in asymptotics. This is a T-reversed state, 
with the corresponding "left" momentum eigenvalue k that is related to that of the "right" eigenstate as k = —k*. 
This agrees with the symmetry properties of the S matrix, 



S(k) = S*(-k*) = S^i-K), 



(41) 
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and assures that left and right energy eigenvalues are complex conjugate. 

Numerically, the decay amplitudes can be calculated directly from (|51Jl or (J3SJ; it has been demonstrated in [3{| for 
proton emitters that these methods are equally effective in practice. The effective non-Hermitian s.p. Hamiltonian can 
be solved resulting in Gamow states via an iterative procedure based on the Green's function, similar to the approach 
discussed in [36|, |32| • Green's function is constructed for the free particle case V = using Coulomb functions at some 
momentum ko and the Siegert boundary conditions, 

g(r,r') = ±F(k r < )O+(k r > ). (42) 
ko 

The r< and r> denote the smaller and the larger of r and r' , respectively. The integral equation for the resonant 
state becomes 

/>OC 

u(r) = / g(r,r') [k 2 - k 2 - 2(iV(r')] u(r')dr'. (43) 
Jo 

This leads to the following equation for the radial part u(r): 

u(r) = if(r) |^°° 0+(r') [k 2 - k 2 - 2[iV] u(r')dr'X + j^O+(r) | jf F(r') [n 2 - k 2 - 2fiV] u(r')dr'^ . (44) 

The complex momentum k is determined self-consistently with the decay flux defined by the outgoing component in 
the above equation. 



B. Threshold behavior 



The behavior of the decay width and the self-energy term A(e) in the vicinity of threshold is particularly important. 
The one-body decay is an instructive example of the CSM at work. Following the above definitions we evaluate 

^■Mrf! = A(e) _i 7(e) (45 ) 

e — €j + lO 2 

in the vicinity of the single-particle threshold at zero energy using Eq. i|34|) which we decompose as 

aj (ej,e) = (j\V\j; e 3 ) + ( £j - e)(j\j; e,). (46) 

Only the first term in (|46|l leads to a pole in Eq. I|45l) . The main contribution comes from low-energy scattering 
states, namely in the limit when the de Broglie wavelength of the scattered particle exceeds the range of the potential. 
For a charged particle, the energy behavior of the amplitude (|34|) follows from that of the regular Coulomb function. 
For a neutral particle, in this limit Ft = (kr) l+1 / (21 + 1)!! oc e^ l+1 ^ 2 . Thus, we can assume that 

(j\V\j;e J )=^(^-y+ 1 / 2 , (47) 



= /So, i mi / r^V{r) U] (r)dr . (48) 



where the constant sej is 

V5F(2I + 1)!I Jo 

In the integral l|45|l only the term oc ee 2 contains the contribution from the pole if e > 0. Direct integration produces 
the result 

A(e) = (i| - Hqq -QV- VQ\j) + e(j\Q\j) + 7rae 2 e(-e) e ^, (49) 

and for imaginary part 

7j -(e) = 2^ae 2 e(e)e' +1 / 2 . (50) 

Here Q = J d€j\jej)(jej\ is the projection operator into Q space; O is the Heaviside step function. The first two 
terms in Eq. (|49|l appear as a correction to energy due to the possible non-orthogonality between the states of Q and 
V. They are small in any reasonably selected situation; they are identically zero in the examples shown below where 
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spaces V and Q are obtained from the full numerical solution of the Woods-Saxon potential. Only the second term 
in l|49|) depends linearly on energy. 

The last term in Eq. Q49J1 appears only below threshold and represents virtual excitation into continuum induced 
by interaction. This contribution is a continuous function of energy, that is also smooth except for the s-wave neutral 
particle decay channel. The behavior of the self energy A(e) including the s-wave ~ y/e type cusp is well studied in 
hadron physics, see also |33l. l39l|. The scaling coefficient eej that enters this expression is typically small; the same 
coefficient determines the energy scaling of the decay width above threshold, Eq. (|50|) . see Fig. EI The behavior 
7 ~ e l + 1 / 2 is consistent with the phase-space volume for one-body decay. It is universal, so that, if the potential is 
being adjusted to reproduce a certain resonance energy, then the behavior of the width as a function of resonance 
energy remains the same [38|, |3{| . 

C. One-body channels in a many-body system 

The one-body decay amplitude in a many-body system is given by the s.p. decay and the spectator overlap, 

At(E)= aj (e) (l;JV|6}|o;JV-l), (51) 

where energy of the continuum state is E = E a + e and c = {a,j}. This amplitude is to be directly used in the 
non-Hermitian effective Hamiltonian (|l(Jfl . Even on the level of s.p. decays, the approach outlined here goes beyond 
the consideration based solely on spectroscopic factors. Here the non-Hermitian part of the effective Hamiltonian is 
not a s.p. operator. Indeed, even s.p. decays can generate significant restructuring inside the nucleus. Effects, such 
as shape changes or changes in pairing coherence, are extremely important for the physics of nuclei far from stability. 
The energy dependence that was discussed earlier is another distinct feature. 

Although in all calculations presented in this work we use a general form of the effective Hamiltonian, below we show 
a set of approximations that establish a correspondence with the traditional SM description of decay. For simplicity 
we assume that all internal s.p. states can be identified by spin and parity, i.e. the space is small enough not to 
include several major shells, so that there is no need in main quantum numbers. The non-Hermitian term in the full 
Hamiltonian coming from s.p. decay channels is then diagonal, 

(i\W(e)\2) = 2k5 12 J2 M^)l 2 l(a;A-%|i ; AOI 2 , (52) 

c (open) 

where s.p. energy satisfies E = 6j + E a ; below we suppress the energy argument if it is unambiguous. 

In the case of remote thresholds, the decay amplitudes become essentially independent of energy, and the set of 
continuum channels c = {a,j} includes almost all possible daughter states a. The completeness in a and energy 
independence can be then used to simplify (|52[1 . 

(1\W\2) = 2 7 n5 12 £> i | 2 (1; N^l; N). (53) 

3 

As a result, W becomes a s.p. operator that assigns a width l|36[l to each s.p. state j coupled to continuum, 
W = 52 . jj bjbj. The same picture emerges when the residual SM interaction is weak. Then s.p. motion masters 
the dynamics, and the sum over daughter systems a in Eq. 152(1 is dominated by a single term. This again leads 
to Eq. (|53[l that is energy dependent, W(E) — ^jlj{E)b^bj. The operator W here can be conveniently combined 
with the SM Hamiltonian just by introducing complex s.p. energies for unstable orbitals in the mean field (a simple 
instructive example was shown in 40] ). Another point to be mentioned here is related to the treatment of the 
non-Hermitian part. In the full CSM diagonalization, virtual transitions to continuum and real decays influence the 
internal structure. This is particularly important at strong continuum coupling when coherence with respect to decay 
leads to the super-radiance phenomenon |7l.l25ll41|. 

The second situation leading to the SM picture is the limit of weak continuum coupling when the matrix W can 
be treated perturbatively. In the lowest order, decays do not affect the internal state and we can solve the Hcrmitian 
problem first to obtain a many-body state \a) with the width given by the expectation value 

T a = {a\W\a)= 1] T ] {a). (54) 

This expresses a many-body decay width as a product of the s.p. width and the spectroscopic factor Tj(l) = 
(l;Ar|6t&.|i ;i v>. 
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D. Two-nucleon emission 



A two-body decay channel state is fixed by an N — 2 nucleus in its eigenstate a and a state of two nucleons in the 
continuum, Eq. Q2H[). These states are characterized by total energy E = E a + €j + 6j> combined of the daughter 
energy E a and energies of emitted nucleons in s.p. channels j and j', total angular momentum and isospin of the 
emitted pair. Unless we are dealing with a bound two-particle continuum state, the channel, besides total energy E, 
has another continuous index describing the energy distribution between the particles. 

As earlier, the two-body transition amplitude is generated by the matrix element (l\H\c; E) . Two different con- 
tributions can be identified as "direct" and "sequential", Fig. 2. The defined above one-body part of the total 
Hamiltonian H cannot contribute to the direct decay vertex 1 — * a + j + f '. The two-body interaction responsible for 
this direct transition is discussed below. Even without two-body interactions, the "dressed" vertex is not zero since 
the N — 2 daughter state is a part of the virtual cloud of the N — 1 system to which one-body transitions are allowed. 
Such a second order perturbative sequential mechanism of decay is possible regardless of whether one-body channel 
is open or closed. Similar a ppr oach and classification of processes have been recently discussed in the context of 
two-proton radioactivity 42, 43] where on one hand the case is complicated by being a three-body Coulomb problem 
[44|. while on the other hand extremely weak decays do not affect the internal nuclear structure allowing the use of 
the traditional real energy SM . 



e: 




FIG. 2: (Color online) Diagrams for two-body decays, sequential (left part) and direct (right part). 



1. Sequential decay 

The dressed s.p. vertex for two-nucleon emission (1\H |c; E) can be calculated with the help of Eq. © and solution 
©, 

„ A\ m (Ea + e, E) A^ aj ' } (E a + e', Eg) , \ 

Assuming that internal and external s.p. states are orthogonal and using one-body decay amplitude (|51|l for the 
second order process that takes place via a virtual state (3 being suppressed by the energy barrier, we obtain 



A c i{E) =}^a j (e)a j ,(e') - +£ 3 (je,j} «-> {e',/}J 



(56) 



where we are considering only on-shell decay so that the initial state with energy E = E a +e + e' undergoes a transition 
via the intermediate state with energy Ep + e. This is reflected in the denominator Sg + e = Ep — E + e, where Sg 
stands for separation energy from the intermediate system. Since c here contains a continuous index, the summation 
over channels is expressed as 



(l\W{E)\2) =n J dede' 5{E-E a -e-e')A L 1 {E)A c 2 *{E). 



(57) 



In order to avoid double counting due to fermion permutation we included a factor of 1/2 while making domain of 
integration symmetric. Eq. Ij57(l may contain poles in open one-body channels corresponding to real states of an 
N — 1 system through which the two-body decay process can proceed. Such processes can be viewed as a second 
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order correction to the one-body decay via the daughter state as a resonant state. Here we discuss only the off-shell 
contribution. An example of sequential decay is shown in Fig. [3 

The calculation of the non-Hermitian term W in the most general case was carried out numerically in the examples 
shown in the Sec. IIVI Here we illustrate the case of a weakly decaying state |otj) when the decay width can be 
approximate d by the diagonal element, T = (cti\W\ai). In this limit the problem is close to that of sequential decay 
discussed in [43, |45j • 

Taking into account the direct and exchange contributions to Eq. I|56|l , assuming spherical symmetry and performing 
the summation over magnetic quantum numbers, we obtain two spectroscopic factors, 

= S »' S ^' (2a i + l)(2/3 + l) ' (58) 

and 

TX _ , u p+p< J on j 0' \ («/IIVII/3 , )*(/?1l& J ll^)*(a/||6 J ||/3)(/3||b / ||a a ) 

T ( 1} \a f f J 2^T1 ' { ' 

where (and below) we use the Greek index of the many-body state as a symbol of angular momentum. The reduced 
matrix elements are defined as in |32T |. The total width in the channel |a,-) — > |a/) is given by 

T-^f^m -«-«') E,*<«h,w + ( 8, + ^ + .) ) <"» 

As an illustration, we discuss a case of the + — > + two- neutron decay. Here T x = T d = T, and we can view a 
full width as a sum of partial widths that depend on the final state a t and intermediate state 0. The characteristic 
energies in the problem are the one-neutron separation energy to the state 0, Sp = Ep — E ai , and two-neutron decay 
energy, E t — E ai — E a . (minus two-neutron separation energy). We assume the situation with small available energy 
E t = e + e', when the energy scaling of the s.p. decay widths, jj = 2n\a j (e)\ 2 ~ e' +1 / 2 , can be used. Introducing 
q = E t / Sf3 we obtain 

r m n\ ^rll{ E t)li'{Et)E t „ , . , . 

r /3 (0^0) = T B w {q), (61) 

which includes the phase space integral 

(9) = — y (i + ^[i +9( i-,)]^ ■ (62) 

In the limit S n ^> E t we obtain Bw(q) = B(l + 3/2, Z' + 3/2) /ir, where B(x, y) is a beta function. The decay rate is 
suppressed by the energy denominator The two-body decay width scales with energy as Ef. This is consistent 

with the phase space volume estimate, T ~ <5 3 (J^ — P)5(E t — E) Yiiid 3 h / £i) , where P is total momentum and the 
product runs over the fragment indices including the daughter nucleus. For the isotropic case integrated over all 
angles, the width is proportional to 7 ~ k ~ E^ 2 in one-body decay; the same assumptions lead to T ~ E? for the 
two-body decay (three-body final phase space). 

In the opposite limit, S n <C E t or q — > 00, we have Bw(q) — B(l — 1/2, Z' + l/2)/(27rq 2 ). This expression diverges 
for an s-wave. The exact integration in the s-wave case gives 

r _ fo -► 0) = il{ E t)~£pE t . . 

4(25^ + Et) y/SpiSp + Et)' 

Note the divergence in the s-wave channel when separation energy goes to zero. The s-wave state in the intermediate 
nucleus is so broad that even being slightly higher in energy it still poses no barrier for the sequential decay. The 
generic behavior of the isotropic decay width as a function of energy can be traced from Eq. JfJ^J), where T ~ Ef 
if energy is low, E t <C Sp, but once energy is getting above Sp the behavior changes to T ~ \ftS~t- This shows that 
the presence of the one-body resonance changes the nature of sequential decay reflecting the one-body phase space 
characteristics. 



12 



2. Direct decay 

Direct two-body transitions are generated by the two-body part of the Hamiltonian Hvq that takes a nucleon pair 
coupled to angular momentum L, pr, = [bjbj/]i,, from internal space V and transfers it to the two-body continuum Q, 
see the right part in Fig [3 The transition amplitude has a generic form, 

A c 1 (E) = a^\e 1 ,e 2 ) (1; N\ (pf^ \a; N - 2) , (64) 

where a direct two-body transition amplitude, a« lJ ' a ^(ei, e 2 ), is introduced (not to be confused with s.p. amplitudes 
a J (e)). This amplitude can be calculated for a given two-body interaction. For example, assuming for simplicity a 
coordinate form V^(r, r'), where r and r' are particle coordinates in the mean-field frame, the amplitude can be 
expressed following the definition of Eq. 0) and normalized free particle states Fj as 

o«^)( Ciea ) = (j 1 j 2 \V^\j 1 ,e 1 ;j 2 ,e 2 ) = — %= f°° drdr' F h {r)F h {r')V^ {r,r')u h {r)u h {r') . (65) 

7rv k\ k 2 J 

The low-energy behavior of the direct amplitude can be understood without specification of the residual interaction 
by taking into consideration the long wavelength behavior of the Bessel functions associated with the regular solution, 
Fi(r) ~ (kr) l+1 at k -> 0. The decay rate can be estimated by integration over continuous channel variables as in 

(53, 

S(Ek - £l - e 2 ) dt x de 2 ~ E l t 1+h+2 . (66) 

The same answer as the one obtained for the sequential transition reflects the nature of three-body final phase space. 
The direct transition however is not suppressed by the one-body energy barrier and is not related to decay amplitudes 
for one-body decays. 

We conclude the discussion of two-body decay processes with the word of caution: the full amplitude in a given 
decay channel is the sum of the direct and sequential contributions and the observed width or cross section carries 
their interference. 

IV. APPLICATIONS 
A. How the method works 

Before turning to the discussion of the specific results we outline the stages involved in the calculation. A special 
feature of our approach is the exact treatment of threshold behavior. This requires the knowledge of the topography 
of thresholds and therefore has to rely on the preceding solution for the daughter systems. In this way we come to 
problems of considering the daughter chains in their entirety. Thus, except for the case of the traditional SM, where 
the Hcrmitian Hamiltonian matrix is diagonalized separately for each nucleus, in all calculations the nuclides (in our 
examples below the isotopes) are coupled by the decay chains. 

The procedure starts from the closed core (specifically, 4 He or 16 0) and continues toward heavier isotopes so that 
the properties of all possible daughter nuclei are known prior to each new calculation. The process of calculating 
resonant states is iterative. We start with a given state obtained from the conventional SM. For this state we review 
all possible decays, including the two-body ones. By conducting scattering calculations at relevant energy or via 
appropriate power-law dependence of Eq. (|50|l at low energies and including the spectroscopic factors from structural 
rearrangement of spectator nucleons we determine the contributions from each decay channel to the non-Hermitian 
term W. The diagonalization of the full Hamiltonian with the non-Hermitian part results in the next iteration for 
energy and width of the state under consideration. For each resonant state this process continues until convergence 
is reached. The Breit-Wigner definition of resonances is used so that scattering calculations are done at real energy. 

By construction of the model, with A(E), Eq. (11), assumed to be included in the adjusted SM Hamiltonian 
with neglected energy dependence, bound levels coincide with those in the standard SM. For unbound states, this 
choice makes our internal propagator equivalent to the i?-matrix used in spectroscopic analysis of experimental data. 
Therefore this design is most suitable for the effective extraction of interaction parameters from experiment. 

The continuum coupling restructures internal states, and energies of resonances above the decay threshold deviate 
from the SM predictions. For narrow states and well separated resonances, the resulting effect is small. The spherical 
shape of the semi-magic system is stable; in addition, strong collective pairing reduces the effect of decays onto internal 
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FIG. 3: Single-particle decay width as a function of resonance energy for the p and sd shell models. The solid lines represent 
results of the reaction calculations based on the Woods-Saxon potential appropriate for the corresponding case (helium or 
oxygen) discussed in text. One experimental point on the p-state corresponds to the ground state 3/2~ resonance in J He 
(e = 0.895 MeV and 7 = 648 keV). Two points for the d-state correspond to experimentally known 3/2 + resonances in lr O 
(e = 0.941 MeV and 7 = 98 ± 5 keV) and 19 (e = 1.540 MeV and 7 = 320 ± 25 keV The near-threshold power law fits 

of Eq. (15011 are shown with dashed lines. 



structure in these particular examples. Whence, for s.p. decays, the case (a) in Table [U the use of the spectroscopic 
factor approximation discussed in Sec. IIII CI is rather good. Here we are far from strong coupling to continuum 
|l2t I25L IHa l54j that may cause an internal phase transition with formation of broad (super-radiant) and very narrow 
(trapped) states. In nuclear physics this phenomenon separates compound and direct reactions [2^, [5^, |5(| . A trace of 
this effect is seen in Table where decaying states have their resonant energy shifted from SM prediction, the lowest 
states of given quantum numbers are pushed into continuum ■ 

The cross section calculations follow directly the formalism outlined in Sec. Ill HI The large scale repetitive matrix 
inversion and instability in the vicinity of the poles presents a significant technical challenge. In this work we 
employed a new numerical method that expands the unperturbed propagator l|16|l in the time domain using Chebyshev 
polynomials. The entire procedure is similar to the Lanczos technique and involves only matrix vector multiplication 
- a fast operation with sparse matrices. The resulting R matrix is then used in the Woodbury equation (|18f) leading 
to a full propagator and scattering cross section. 



B. Helium isotopes 

For the chain of helium isotopes from 4 He to 10 He, the results are summarized in Table. [I] The internal valence 
■p-space contains two s.p. levels, P3/2 and P1/2, the a-particle core is kept inert. The effective interaction within 
this model space and s.p. energies are borrowed from |46L l47j. Without additional terms, this would be merely a 
conventional SM leading to the bound states that are listed as E(SM) in Table [I] 

The one-body part of the coupling Hamiltonian H-pq is defined using the Woods-Saxon potential. However, at 
low energies the s.p. decay widths can be well approximated by near-threshold dependence of Eq. (41) as 73/2(6) = 
1.08 e 3 / 2 MeV and 7172(e) = 0.151 e 3 / 2 MeV for p 3 / 2 an d P1/2 states, respectively, see Fig. [3] The results of calculation 
limited to s.p. decays are shown as case (a) in Table [I] The "Borromean" 6 He nucleus requires two-body decays. 
The sequential two-body decay as a second order process built on the one-body amplitudes involves no additional 
parameters. This process is included in the results shown as case (b) in Table [i] 
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3/2 


0.992 


0.992 


0.992 


0.992 


0.895 


0.6 


0.6 


0.6 


0.648 


5 


1/2 


4.932 


4.932 


4.932 


4.932 


4.895 


4 


4 


4 


4.1 


6 





-1.379 


-1.379 


-1.379 


-1.379 


-0.973 














6 


2 


0.515 


0.515 


0.529 


0.529 


0.825 





0.248 


0.248 


0.113 


6 


2 


4.745 


5.25 


5.25 


5.25 




2.566 


2.566 


2.566 




6 


1 


5.889 


5.32 


5.32 


5.32 




0.922 


0.922 


0.922 




6 





11.088 


10.911 


10.911 


10.803 




5.532 


5.532 


12.303 


broad 


7 


3/2 


-1.016 


-1.016 


-1.016 


-1.016 


-0.528 


0.046 


0.046 


0.046 


0.15 


7 


1/2 


2.24 


2.239 


2.253 


2.253 


2.172 


2.357 


2.689 


2.69 


3 


7 


5/2 


2.85 


2.888 


2.911 


2.911 


2.393 


0.727 


0.944 


0.944 


1.99 


7 


3/2 


4.495 


4.379 


4.222 


4.22 


(5.273) 


0.541 


1.113 


1.246 
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7 


3/2 


10.223 


8.857 
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7.818 


16.379 
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-3.591 


-3.591 


-3.591 


-3.591 
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8 


2 


0.19 


0.196 


0.191 


0.19 


-0.308 


0.231 


0.506 


0.53 




8 


1 


2.427 


2.304 


2.331 


2.321 




1.026 


1.455 


1.418 




8 





6.376 


6.003 


6.527 


6.489 




5.286 


3.456 


15.449 
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2 


6.882 


6.839 


6.538 


6.572 




2.283 


13.86 


14.94 
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1/2 


-1.992 


-1.992 


-1.992 


-1.992 


-1.958 


0.634 


0.634 


0.634 


0.1 


9 


3/2 


2.805 


2.801 


2.802 


2.797 


-0.6 


1.557 


2.425 


2.443 


0.8 


10 





-1.649 


-1.649 


-1.649 


-1.649 




0.073 


0.504 


0.746 


0.3 



TABLE I: Comparison of conventional SM and CSM with data for He isotopes (all numbers in MeV; energies are measured 
from the ground state of 4 He). The first two columns indicate the mass number and spin of the state. The next five columns 
compare energies as follows: _E(SM) - traditional shell model; E(a) - CSM with only one-body decays included; E(b) - CSM 
with one-body decay and its second order contribution to the two-body process; U(CSM) - full CSM including the direct two- 
body decay mode; _E(EX) - experimental data (some of them have large uncertainties and depend on the way of analysis). Last 
columns compare decay widths from CSM calculations with data |50ll5lLl52j|: the SM calculation gives only discrete energies. 

The judgment about the direct two-body decays, their importance, and whether two-body terms have to be included 
in H-pQ can be made with comparison of theory versus experimental data. Considering the simplest two-body case 
of 6 He we observe that the included sequential decay provides a reasonable description for the lowest 2^ resonance in 
6 He. However, the need to include direct two-body decay comes from consideration of a broad asymmetric resonance 
structure at excitation energy ~ 3 — 15 MeV |4Sll4l^| ascribed to neutron skin oscillations. Although an experimental 
controversy exists in relation to this structure, our results below, as well as earlier studies, see |48|. suggest a direct 
two-body breakup process. Based on the experimental exclusion of 2^~ possibility by |49| , we tentatively identify this 
resonance as a spin zero pair-vibration excitation generated by a coherent transition of the L = neutron pair between 
the two s.p. orbitals and introduce direct two-body decay with the pair emission in the L — channel. At this level 
of experimental precision and with the sequential decays being already included, there is no compelling evidence for 
direct decays in other angular momentum channels. 

Without introducing a specific isotropic form of the two-body interaction for H-pq we adopt a low-energy "pairing" 
approximation assuming that all L = neutron pairs couple to continuum with the same amplitude a^ i=0 ^(ei, t 2 ) = 
(ei + e 2 )/3v / 27r. The numerical constant here sets the strength of the residual two-body interaction. This is the only 
parameter of the model, and with limited experimental knowledge it was selected to roughly reproduce the broad 
resonant structure in 6 He. The results of the full calculation are shown in Fig. 0] and listed as (CSM) in Table [I] 
Fig. 0] also shows elastic neutron scattering cross section from the ground state of the N — 1 nucleus in M = (even) 
or M = 1/2 (odd mass) magnetic moment state. The cross section is computed with the same CSM parameters but 
within the case (a) discussed above. 

Below we comment on some of the features of the results. 

(i) The ground states of 4:6,8 He are nucleon-stable in agreement with experiment. 

(ii) The results reveal information about structure and dominant decay modes. For example, for the 7 He isotope, 
that has not been used for an adjustment of interactions, they agree with recent experiments |5ll l52|. Our results 
support the "unusual structure" of the 5/2 ~ state identified by 52]. Due to its relatively high spin, this state, unlike 
the neighboring l/2~ state, decays mainly to the 2 + excited state in 6 He. Concerning the deviations from the data 
in the width of this state and 3/2~ state, we may comment that experimentally although there are still serious 
uncertainties, it is most likely that the contribution from the triple decay sequence, sequential two-body followed by 
one-body, is important for the lifetime of the 3/2~ state. The processes of the third order one-body decay or second 
order combination of direct two-body followed by one-body, are not included in this calculation but in light of the 
experimental data they should be considered in the future. 
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FIG. 4: (Color online) CSM results for He isotopes. The states in the chain of isotopes starting from 4 He (top) to 10 He (bottom) 
are shown as a function of the energy relative to 4 He. The horizontal dotted lines separate each isotope and for each case states 
from CSM are shown above experimentally observed states. The decay width (in MeV) is shown for each state along with spin 
and parity. The solid lines above CSM states shows the elastic neutron scattering cross section from the spin polarized state of 
TV — 1 isotope, see text for details. 



(iii) Our results and, most importantly, the theoretical interpretation of the role that one-body and two-body 
interaction terms play in the dynamics of decay agree with earlier findings in |Tflj | where a different method involving 
a full diagonalization without separation of continuum into channels was used. 

(iv) For the heaviest isotopes and, in particular, for the discussion of the presently interesting case of 9 He, the 
extension to the sd shell is necessary. 

(vii) The scattering cross section is non-zero only above the thresholds (the ground state of the A — 1 nucleus). 
Cross section curves are generally not symmetric and have neither Gaussian nor Lorentzian shape. For low-lying 
states with the width big enough to reach threshold, the distortion is particularly noticeable. 

(viii) For broad resonances, the identification of resonances (or poles of the scattering matrix) with a peak in the 
cross section is ambiguous. Interference between different resonances, including rather remote once, is significant in 
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"the Q- value is measured for the excited state in daughter system. 

TABLE II: Lowest resonance states in the chain of oxygen isotopes. The experimental data on the left, (EXP) - energy of the 
state (MeV), Q - energy above threshold (MeV), V - width (keV), are compared to the theoretical results on the right. The 
decay mode in the second column indicates the decay branches assumed by experimentalists [Hoi Ir3(i|| . 

this case. 



C. Oxygen isotopes 

The internal space here is represented by the S1/2, d 3 / 2 and d 5 / 2 s.p. orbitals composing the usual sd-shell model. 
The standard SM interactions (USD 30] or its slightly modified version for heavier isotopes HBUSD were used 
in these calculations. 

Since we use the Woods-Saxon potential for reaction calculations, it is important to make sure that the potential 
create proper resonant states. For few low-lying s.p. resonances with large spectroscopic factors, we adjust the 
depth of the potential so that the correct s.p. resonance energy is indeed ensured. The Woods-Saxon potential 
parameterization with the mass dependent depth j58L 159| was demonstrated to reproduce the s.p. resonances and 
bound states with good precision for nuclei around mass A = 16. For the majority of high-lying states, the potential 
is not readjusted, here the decay amplitude is computed directly from Eq. I|34() . The resonance phenomenon in this 
case is due to the many-body effect where the complexity of the many-body states makes the overlap in Eq. I|51|) 
small. 

In Fig. [21 the widths of the resonant states S1/2, d 3 / 2 and 1^5/2 are shown as a function of their energy being found 
with the aid of the Woods-Saxon potential with variable depth. The curves are limited to the near-threshold region 
approximately determined as kR < I (and kR <C 1 for I = 0) where R is the nuclear radius. For oxygen, this limits 
the d-wave at about 3 MeV. The curves are close to straight lines displaying the appropriate power law scaling 15UH of 
the decay width as a function of energy. The lines can be fit by equations (with e in units of MeV) 7(51/2) = 16e 1//2 , 
7(d 3/2 ) = 0.15e 5 / 2 , and j(d 5/2 ) = 0.04e 5 / 2 . These fits are shown with dashed lines. 

In Fig. [3] we show an overview of the full CSM calculation for oxygen isotopes within the sd shell. Some of the first 
results were reported earlier 0, HI ■ 

Our calculations are in a good agreement with available experimental data for the oxygen isotopes |50t l6fl | . One 
has to emphasize again that, with the assumption of the self energy term A being a part of the conventional SM 
Hamiltonian, the bound states are exactly the same as obtained in the usual SM calculation. The novelty appears 
when the states above decay thresholds are considered. The properties of a few experimentally identified resonances 
in oxygen isotopes are compared to the CSM predictions in Table ITD 

There are two cases in 17 and 19 corresponding to the neutron emission from the d 3 / 2 orbital. These rather 
pure s.p. decays involve little many-body physics and are well described with the Woods-Saxon potential model, Fig. 
13 Here the SM may be used to predict the neutron decay energy, however there is an obvious improvement if the 
experimental Q- value is used in the s.p. reaction calculations. The power law for the energy dependence of the decay 
width can enhance the predictive power of the description. 

The two-body case of 18 is more complicated, however there are only 14 states in the sd shell model. The 
comparison of the level structure and neutron scattering cross section with data is shown in Fig. [5] The SM below 
neutron decay threshold (zero on the plot) is in an impressive agreement with the data. Here, for all sd SM states 
shown with longer ticks (in red, color online), the experimental counterpart can be easily identified in the observed 
spectrum. Practically all remaining observed negative parity states below neutron threshold can be identified with 
particle-hole excitations using an extended p — sd — pf SM. In theoretical calculation, these states are shown with 
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FIG. 5: (Color online) CSM calculation for oxygen isotopes with the HBUSD interaction. The states (vertical bars) are plotted 
as a function of absolute energy computed relative to ground state 16 O. States shown with solid black lines are stable in our 
model; they are either below decay thresholds or with decays forbidden due to the angular momentum restrictions in the 
selected valence space. States from yellow (color online)/lighter shade of gray (long lifetime) to red/darker gray (short lifetime) 
are resonance states. For some of the low-lying states the energy (in units of keV), spin and parity are given. Along with the 
resonant level structure above thresholds, the elastic d-wave neutron scattering cross section, off a ground state of a daughter 
nucleus with M = (even mass) or M = 1/2 (odd mass) magnetic quantum number, is plotted as a function of energy. The 
cross section is obtained as a part of the same CSM calculation. 
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FIG. 6: (Color online) The lower panel is the result of the CSM calculation of the level structure and neutron scattering cross 
section from the ground state of ir O. Marked levels corresponding to the sd-shell model states are included in CSM calculation. 
All other levels result from the expanded sd — pf shell calculation (see text); the configurations of this type are not coupled 
to continuum and corresponding (negative parity) states do not appear in the cross section plot. The upper panel corresponds 
to the empirical level scheme in ls O. The neutron scattering cross section on the upper panel is indirectly inferred from the 
14 C(an) 17 reaction |61f and other compilations 60, 62. 63]. The parts of the cross section shown with solid line are expected 
to be dominated by the positive parity / = or I = 2 partial waves and thus can be distantly juxtaposed to the curve on the 
lower plot. 

small lines (the WBP cross-shell interaction from Ref. [3]| was used). 

Due to the resonant nature of states, the experimental picture above neutron decay threshold becomes more 
ambiguous. The lowest near-threshold state 2 + may correspond to one of the higher-lying 2 + states in our model, sec 
Fig. [fJland Table ITU Due to the significant difference in energy above threshold between this state and its possible 
theoretical counterpart, comparison of the decay widths is inappropriate. Following this state, there are rather narrow 
4+ and 1 + states (the experimental assignment of 1 + for the latter state is still uncertain), see Table ITT1 and Fig. 
Although in theoretical calculation these states appear in a different order their observed and calculated widths appear 
to agree. Beyond this point there are several broad 100-200 keV unidentified resonances observed in experiment. They 
may be juxtaposed to 150-400 keV wide 2 + , 3 + , and 1 + states appearing in the CSM calculation. To support this 
argument, we show on the lower panel of Fig. Hjjthe CSM neutron scattering cross section computed using Eq. ifHJl. 
Since this calculation takes into account only I = and I = 2 partial waves, the peaks from negative parity states do 
not ap pea r. On the upper plot we roughly reproduce the experimental cross section obtained from indirect reactions, 
Refs. [(MJ, l6ll I62L l63l | . The portion of the curve corresponding to positive parity resonances that are expected to be 
present in the theoretical calculation is shown by the solid black curve. There is a reasonable similarity in the resonant 
structure corresponding to the set of wide states 2 + , 3 + , and 1 + . In the experiment, the spin-parity assignments for 
the states shown in Fig. |B|in brackets have not yet been confirmed. 

The cross section picture indicates that our interpretation is plausible, while differences, by the factor of 2 to 4, in 
the resonance parameters may be due to different definitions of the width used in the data analysis. The broadness of 
these resonances causing difficulties in experimental analysis and interpretation of peaks is here a trivial consequence 
of the simple two-particle structure of 18 that implies high spectroscopic factors for s.p. decay. The situation changes 
in the next isotope 19 where there is only one broad s.p. 3/2 + state that stands out in the calculations of the cross 
section. The majority of other states are narrow due to the many-body complexity. They would not be visible on 
the plot of the cross section; on the other hand, this allows for the direct comparison of widths and energies shown 
in Table [H] The level scheme inferred from experiments with 19 O fits well the overall picture of oxygen isotopes as 
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shown in the insert of Fig. 

Only little is known about the heavier oxygen isotopes although the situation is going to change in the near future 
with the new radioactive beam experiments. In Table ITT1 we quote some, in our view the most interesting, predictions 
of the model. The cases of 26 and 28 are just beyond the drip line. The ground states of these nuclei are unstable 
with respect to the neutron pair emission. Sec. IIII Dl while being centered around interactions, discussed how to 
tackle such cases. Unfortunately, the uncertainty in the effective interactions complicates the job of predicting. We 
have no firm knowledge about the interaction coupling internal states to the continuum and, unlike the helium case, 
in the oxygen data we could not find a case to determine this coupling phenomenologically. Thus we only consider 
here the case of sequential decay driven by what we believe to be a well defined one-body potential. 

The second problem is that the traditional SM adjusted to experimental data near the stability line cannot be 
extrapolated with full certainty to the vicinity of drip lines. The role of the self-energy term A(E) is another 
question. Indeed, the well established USD interaction predicts 26 O to be bound which is known from experiments 
not to be the case. In our calculations we had to resort to the HBUSD interaction specifically adjusted to heavier 
isotopes. The Q- va hie coming from this interaction in the case of 26 is only 21 keV while the typical SM uncertainty 
in level energies is about 200 keV. This uncertainty entering the power law scaling of decay width versus energy makes 
the lifetime predictions unreliable. 

Thus, while collecting our calculations in Table ITT1 we see our strongest prediction through Eq. (|63|l . In both cases 
of 24 and 26 0, the lowest 3/2+, 1/2+ and 5/2+ states in the adjacent isotopes 23 and 25 are the main candidates 
for the intermediate states involved in the sequential two-body decay. This is because of the low energy barrier in 
combination with large spectroscopic factors of these mainly s.p. states. At the same time, in the case of 24 we 
have separation energies S 3 / 2 = 0.98, S]/ 2 = 4.39 and S 5 / 2 = 5.54 MeV; the corresponding spectroscopic factors are 
T 3 / 2 = 0.28, T1/2 = 0.019 and T 5 / 2 = 2.2 x 10~ 4 . In addition to spectroscopic information, the kinematics of the 
phase-space volume is another essential factor. A simple estimate with the aid of Eq. (|63[) shows that the transition 
through the s-state would dominate throughout the entire region below s.p. threshold. Due to level energetics, two- 
body decay remains significant even in the presence of the open one-body channel. However, the simple power law 
scaling used to obtain the phase space integrals Eq. I|t)2|) may no longer be valid high above threshold. If the ground 
state in 24 were at E t = 1 MeV above two-neutron decay threshold, while still right at the opening of one-neutron 
decay channel into the 3/2+ state of 23 0, S3/2 = 0, the decay width coming from sequential decay would be about 
30 keV, as opposed to 20 eV in Table HJ quoted for E t = 0.021 MeV. 

V. CONCLUSION 

The first goal of this work is to present a systematic and detailed discussion of the continuum shell model as a 
step in the direction of unifying the nuclear structure with nuclear reactions. We amplify and extend the ideas and 
methods started in our earlier works 0, @ ■ In this presentation we clarify the CSM formalism and show its relation to 
the standard SM. One of the important points of this work was physical interpretation of various results. On the shell 
model side, we highlighted the meaning of solutions of the energy-dependent non-Hermitian effective Hamiltonian and 
identified the procedures to be taken in relation to different definitions of resonant states. On the reaction side, we 
show how the scattering matrix, cross sections and related quantities can be calculated; the unitarity properties of 
the scattering matrix built in the model are emphasized. 

The general discussion of s.p. decays centers around the potential s.p. problem. This textbook problem is not only 
a central part of calculations but it also provides an important parallel to the full CSM description. Here the Gamow 
states, decay amplitudes and scattering matrix can be calculated through numerical solutions of the Schrodinger 
equation in coordinate space. This allows one to establish a transparent relation and interpretation of the same 
quantities in the full CSM. 

The consideration of two-body decays is one of the significant advances of the present formulation. Considering 
the one- and two-body terms in the part of the Hamiltonian that links the internal shell model and external reaction 
space we get keys to the sequential and direct decays, respectively. In the case of sequential decay, a second order 
one-body process, we discuss the transition through the resonance tail, namely the role of broad one-body resonances 
located in the intermediate nucleus above threshold. The direct, or correlated, decay processes are strongly related 
to the problem of pairing and other coherent effects in the continuum. This problem, important also for physics of 
neutron stars, is still far from being solved. 

The last section of this work shows practical applications of the model. The self-consistency in energy, proper 
open channels and realistic reaction calculations, parent-daughter structure relations through the decay chain, — arc 
discussed as the essential elements of the model. The diagonalization of the full Hamiltonian with both Hermitian and 
non-Hermitian parts is emphasized as an important component for treating properly the mutual influence of structure 
and reactions. The discussion of the extreme effects of this nature, such as super-radiance, is outside the scope of this 
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paper, see pi IHl Ip] . 

Another goal of this work is to report the advance in practical applications of the CSM and in particular to 
demonstrate new cross section calculations performed in the same unified framework and presented along with the 
bound states and resonance parameterizations. Comparison with experimental data indicates a satisfactory agreement 
for both helium and oxygen isotopes. The possible role of sequential decay in heavy oxygen isotopes is discussed and 
predictions for the decay width are given. Coming experimental data will be instrumental for further development of 
theory. 

We did not discuss here the computational problems which are of a higher level of difficulty compared to the normal 
SM. In fact, the novel methods for constructing Green's functions for large-scale calculations using the Chebyshev 
polynomial expansion and Woodbury equation stand behind the presented results. However, the technical and nu- 
merical details are outside the scope of this article; they will be presented elsewhere. From the conceptual viewpoint, 
the basic question of effective interactions remains unsolved. We used here a semi-empirical method of combining 
the SM experience with finding the missing cross-space matrix elements from the solution of the scattering problem 
and numerical fits based on general requirements of quantum-mechanical threshold behavior. The inclusion of giant 
resonances, more complicated decay modes and cluster channels is also on the agenda for future work. 
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